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Nonvanishing anisotropic magnetoresistance in Rashba two- 
dimensional electron systems with nonmagnetic disorders 
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PACS 73 . 43 . Qt - Magnetoresistance 

PACS 72 . 25 . Dc - Spin polarized transport in semiconductors 
PACS 71.70.Ej - Spin-orbit coupling 

Abstract. - We study anisotropic magnetoresistance (AMR) in a spin-polarized two-dimensional 
electron gas with Rashba spin-orbit coupling and nonmagnetic disorder collision. We show that 
AMR exists, arising from the combined effect of in-plane magnetization, spin-orbit coupling, and 
nonmagnetic remote disorder scattering. Further, numerical evaluation demonstrates that the 
smoothness of the remote disorder can strongly affect AMR, and this AMR is sensitive to the 
electron density. Large magnitude of AMR (« 24%) is obtained for low density system with 
strong spin-orbit splitting. 
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Introduction. — Magnetotransport phenomena in 
ferromagnetic semiconductor, such as anisotropic magne- 
toresistance (AMR) [1-12], have attracted significant at- 
tentions due to applications in the emerging field of spin- 
tronics [13,14]. AMRs, including longitudinal and trans- 
verse AMRs, are the response of magnetoresistance to 
the relative angle between magnetization and current in 
magnetic materials. Both the longitudinal and transverse 
conductivities show the symmetric feature: a xx (Mo) = 
q xx (-M ) and a yx (M ) = a yx (~M ), where the magne- 
tization Mo is usually in the two-dimensional plane (x- 
y plane). However, one should note that for ordinary 
charge Hall effect (including anomalous Hall effect [15]) 
the transverse conductivity obeys the antisymmetric rela- 
tion, a yx (Mo) = —o- yx (—Mo). Here Afo is normal to the 
plane. 

Experimentally, AMR has been extensively studied in 
diluted magnetic semiconductors recently. Rushforth et 
al. investigated the physical origin of the noncrystalline 
and crystalline components of AMR in diluted magnetic 
semiconductors [6]. Shin et al. explored the temper- 
ature dependence of AMR in ferromagnetic (Ga,Mn)As 
films [16]. A giant transverse AMR was also observed 
in this ternary ferromagnetic semiconductor (Ga,Mn)As 
[9, 10]. In contrast to the extensive experimental studies 
of AMR, the theoretical interpretation of AMR is rela- 
tively poor. The experimental analyses are usually based 
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on a phenomenological treatment [2,3]. The full Boltz- 
mann theory simulations have been made to study the 
origin of the sources of AMR in p-type magnetic semi- 
conductor [6,8,17]. Kato. et al. analyzed the intrinsic 
AMR in spin-polarized two-dimensional gas (2DEG) with 
Rashba spin-orbit coupling (SOC) [7]. They showed that 
AMR vanishes unless the relaxation time is spin-related. 
Recently, Trushin et al. studied AMR for Rashba or Dres- 
selhaus spin-orbit splitting electron system with polarized 
magnetic impurities [18]. In the above theoretical studies, 
the microscopic mechanism of AMR is considered as due 
to the anisotropic carriers lifetime, arising from the com- 
bined effect of the SOC and the scattering by the polar- 
ized magnetic impurities. The spin-dependent scattering 
is the essential factor in AMR [6-8, 12, 18, 19]. In most 
studies, the nonmagnetic disorder potential is taken as 6- 
form type. However in realistic heterostructure, the elec- 
tron density is not large enough to screen the nonmag- 
netic impurities, where the interaction between electron 
and disorder is long-ranged. Hence, the effect of electron- 
impurity scattering on AMR is far from being understood, 
completely. 

In this paper, we employ the kinetic equation approach 
to investigate AMR in two-dimensional electron system 
in the presence of Rashba-type spin-orbit interaction and 
in-plane magnetization. We show that the combined ef- 
fect of SOC, in-plane magnetization, and nonmagnetic 
long-range impurity can lead to AMR. At the same time, 
numerical evaluation demonstrates the disorder-distance- 
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and electron-density-related feature of AMR. The present 
study may provide another mechanism of AMR in spin- 
orbit interaction 2DEG with in-plane magnetization. The 
nonmagnetic wave-vector-dependent disorders, coupling 
to the SOC and magnetization, induces the magnetotrans- 
port anisotropy. 

Theoretical approach. — We consider a 2DEG con- 
fined in a [001]-grown III-V semiconductor heterostructure 
with Rashba SOC and a homogenous in-plane magnetiza- 
tion Mq. The x and y axes are taken along [100] and 
[010] direction, respectively. Hence, the noninteracting 
one-particle Hamiltonian can be written as 

k 2 

H = h a(z x er) • k - M x a x - AI v a v . (1) 

2m y " 

Here m is the electron effective mass, <x = (a x , a yi <j z ) 
are the Pauli matrices, k = (kcos8k,ksm8k) is the two- 
dimensional electron wave vector, a is the Rashba SOC 
parameter, and M = (M Xl M y ) = M(cos£,sin£) = 
gpsM = <7/iaAfo(cos£, sin£) with g as the effective g- 
factor, as the Bohr magneton, and £ as the angle be- 
tween the magnetization M and [100]-axis. 

The above Hamiltonian (TTJ) can be diagonalized into 
H = U^HUk = diag[ei(fc), £2{k)] in the hclicity basis with 
the help of the following local unitary transformation 
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Here the energy dispersion e fJl (k) 
with £ RM (fc) = yja 2 k 2 + M 2 + 2akM sin(£ 
1, 2 as the helix band index, and 



h + (-l)^ RM(fc) 

M = 



Xk = tan" 



_ x ak sin 8k — M cos£ 
ak cos 8k + M sin £ 



(3) 



Now wc consider the quasi-two-dimensional system is 
driven by a weak dc electric field E along [100] direc- 
tion. Obviously, in order to carry out the evaluation of 
the AMR, it is necessary to determine the matrix elec- 
tron distribution function. The kinetic equation for the 
2x2 matrix distribution function p(k) in the station- 
ary linear response regime can be derived, where the 
elastic electron-impurity scattering is taken into account 
in the self-consistent Born approximation [20-22]. Fol- 
lowing the procedure of these papers, the distribution 
function can be obtained as, p(k) — p^(k) + p^(k) + 
p( 2 \k) 7 with equilibrium distribution function p'°'(fc) = 
diag{riF[ei(fe)], Tip [£2 (&)]}, and np(x) as the Fermi-Dirac 
function. Here p^{k) and p^ 2 \k) are collision-unrelated 
and collision-related matrix distribution functions in the 
first order of electric field, respectively. The collision- 
unrelated distribution function p^lk) is off-diagonal ma- 
trix with the elements given by 



eE d\k 
4erm dk x 



{n F [ei(fc)] - n F [e 2 (k)}}, 



where Eq is the strength of the electric field. This distri- 
bution function is associated with the interband transition 
between two spin-orbit-coupled bands, making no contri- 
bution to charge conductivity. However, it is important for 
spin Hall effect, which results in the collision- independent 
intrinsic spin Hall effect [22-24]. The collision-related dis- 
tribution function p^ (k) is determined by the coupled 
equations 



eE, 



dn v (e k i_i) 
dk x 



7r]T|V(fc-g)| 2 fi wl , 

qfj,' 

[P^(k)-p%(fl)]6[e^k)-^(q)], 

(5) 
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P%(k)-p%(q)]s[e,(k)-s,,( q )]. 

(6) 



Here O w < = 1 + (-1)*'+^ cos(xfc - Xq) and = 

(— l) p sin(xfc — Xq)- V{k — q) is the nonmagnetic impurity 

(2) 

scattering potential. Kcp 12 (k) represents the real part of 
the off-diagonal distribution function p\ 2 (k). One should 
note that here the weak scattering limit is assumed, where 
we restrict ourselves to the leading order of the impurity 

concentration. In this case, the imaginary part of the 

( 2) 

off-diagonal distribution function p 12 (k) can be ignored 
completely [20]. In the above kinetic equations, both the 
interband and the intraband transitions are considered. 

In order to study AMR, it is necessary to evaluate the 
drift velocity. In spin basis, the two in-plane matrix ve- 
locity operators read 



-ia ^ 
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(7) 



(8) 



It is clear that the velocity operators in spin basis are 
independent of the magnetization. Moreover, the expres- 
sions of the velocity operators are the same as the ones 
of semiconductor heterostructure with Rashba spin-orbit 
interaction in the absence of magnetization. However, in 
the hclicity basis, the single-particle operators of velocity 
ii = U^ViUk (i — x, y), rely on the magnetization through 
the energy spectrum and the angle Xk, and are given by 



aei(fe) 

—iasinxk 

dk y 

ia cos Xk 



ia sin Xk 

de 2 (k) 



-ia cos Xk 

de 2 (k) 
dk v 



(9) 
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One find that the off-diagonal elements of in-plane ve- 
locity operators are also nonvanishing in hclicity basis. 
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The corresponding macroscopical drift velocities are ob- 
tained by taking the statistical average over them, t'i = 
-k Tr[p(k)vi], and expressed as 



1 y> gggW (a) 

TV 2—* ^ w 



dki 



(fe). 



(11) 



Here A" is the electron density. It can be seen that the av- 
erage velocities only depend on the diagonal clement of ve- 
locity operators. One should emphasize that in clean limit 
approximation, the imaginary part of the off-diagonal el- 
ement of collision-related distribution function vanishes. 
Hence, the drift velocities only relate to the diagonal ele- 
ments of velocity operators and the diagonal elements of 
distribution function. And the expressions of average ve- 
locities become the same as the usual form of two band 
system without interband coupling. We only need Eq. ([5]) 
to determine the diagonal elements of distribution func- 
tion. The real part of off-diagonal elements, Kep 12 (k), 
is essential for calculation of spin Hall effect [21,22] and 
anomalous Hall effect [20]. The longitudinal and trans- 
verse conductivities arc defined by a xx = Nev x /E$ and 
Cyx = Nevy/Eo, respectively. 

For AMR, one can find that the longitudinal and 
transverse conductivities obey the symmetric relations: 
a xx (M ) = a xx (-M ) and a yx (M Q ) = a yx (-M ). We 
can understand these properties as follows: The eigenen- 
ergy e Al (M , k) and angle Xk(M ) satisfy e AI (-M , -fe) = 



TT + 



Xk(M ) 

S 2 h 



Thcrc- 



Sfi (M ,k) and X -k(~M ) 

fore distribution function satisfies p^(— Mq, — k) 

(2) 

—p^fi^Mo, k). Note that for brevity, the argument Mq for 
eigenenergy, distribution function and so on, is dropped 
elsewhere. When Mq — > — Mq, we make transformation 
fe — ► — k in Eq. (fTTj) . This transformation will not change 
the total integral, hence the conductivities satisfy the sym- 
metric property. This property is in vivid contrast to the 
one of anomalous Hall effect, where the transverse conduc- 
tivity obeys the antisymmetric relation. Now we remark 
these relations from the point of view of the time rever- 
sal symmetry. For AMR, the conductivities are related to 
the diagonal elements of distribution function, which are 
proportional to a momentum-dependent effective trans- 
port relaxation time r tr (see footnotcQ). However, the 
anomalous Hall conductivity relies on the off-diagonal el- 
ements. The off-diagonal elements do not depend on this 
effective transport relaxation time, directly, and the dis- 
order plays only an intermediate role [20]. Under time 
reversal, for AMR, <r xx (-M , -r tr ) = -cr xx (M Q , r tr ) and 
a yx (-M , -T tr ) = -a yx (M , r tr ). By considering the re- 
lations between the conductivities and the effective relax- 
ation time, one can obtain the symmetric feature. While 
for anomalous Hall conductivity a yx (— Mq) — —a yx (Mo) 



1 Actually, the diagonal elements of distribution function are pro- 
portional to the quantity with dimension of time, relating to several 
band-dependent relaxation times [20]. We call this quantity as ef- 
fective transport relaxation time. 
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Fig. 1: The relative fractional change in the longitudinal (a) 
and transverse (b) conductivities as functions of the angle £ 
between magnetization in the plane and [100]-axis for various 
remote impurity distances. The thin wine lines are obtained 
for r5-form short-range electron-disorder collision. The electron 
density N = 1.0 x 10 1 cm - . Rashba constant a = 3.0 x 
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eVm and the magnetization M = 2meV. 



under time reversal, and the antisymmetric relation is ob- 
tained directly. 

Numerical results. — We numerically investigate 
the combined effect of Rashba SOC, magnetization, and 
long-range nonmagnetic electron-impurity scattering on 
the AMR in InAs/InSb heterostructure. The long-range 
electron-impurity collision is considered as due to remote 
nonmagnetic charged impurities separated at a distance s 
from the interface. The potential takes the form V(q) = 

U(q)/ K (q) with \U(q)\ 2 = m ( ^) V 2 ^) 2 [25,26]. 

Here I(q) is the form factor and n(q) is the factor related 
to the Coulomb screening, the expressions of which can 
be found in Rcf. [26]. rij is the density of nonmagnetic 
remote impurity. In the calculation the electron effective 
mass is taken asm = 0.04m o with m as the free electron 
mass. The dielectric constant of InAs k = 15.15. In the 
numerical analysis, we have assumed that £f > M, i.e. 
both the minority and the majority bands are occupied, 
with £p as the Fermi energy. 

We first consider the relative fractional changes in the 
conductivity (RFC), which arc defined by 



RFCx = 
RFCy = 



Act, 



Act 



A 



Act 



where Act t 



^avj AfTya; CT^ 



(12) 
(13) 



a av , and 
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Fig. 3: Dependencies of AMR on SOC constant a for different 
electron densities. The magnetization M — 1.5 meV and the 
remote impurity distance s = 50 nm. The inset shows AMR as 
functions of electron density when Rashba spin-orbit splitting 
a = 4.0 x 10 -11 eVm. 



Fig. 2: AMR is shown as functions of spin-orbit interaction 
parameter for various magnetizations (a) and as functions of 
magnetization for various Rashba coupling parameters (b) at 
fixed distance s = 30 nm. Here the electron concentration N = 
1.0 x 10 11 cm -2 . The thin wine lines indicate the corresponding 
AMR for short-range electron-impurity scattering. 



is the average value of the longitudinal conductivity as 
the magnetization is rotated through 360° with respect to 
[100]-axis. The subscript "min" means the corresponding 
minimum value of the fractional change. It is seen that 
RFC is independent of impurity density. 

In Fig. [TJ the relative fractional changes in the longi- 
tudinal and transverse conductivities are shown as func- 
tions of the angle between magnetization in the plane and 
[100]-axis for this nonmagnetic remote disorder. The cor- 
responding thin wine solid line is obtained for S shape 
short-range electron-disorder collision. V(q) — Vo, inde- 
pendent of momentum. It is clear that the longitudinal 
conductivity shows the strong anisotropy for magnetiza- 
tion aligned along various direction when the nonmagnetic 
disorder is long-ranged. However, the anisotropy com- 
pletely vanishes for short-range electron-impurity scatter- 
ing, in agreement with the previous studies [7,18]. The de- 
gree of this anisotropy depends strongly on the smoothness 
of the remote disorder. With the rise of the impurity dis- 
tance, the degree of the anisotropy first enhances, and then 
drops when the distance is large enough. In Fig. Q] (b), it 
is seen that the transverse conductivity also indicates the 
anisotropy for various direction of magnetization. The re- 
mote disorders can affect the degree of the anisotropy of 
transverse conductivity, similar to the longitudinal con- 
ductivity. It is also found that, when the nonmagnetic 



disorder becomes short-ranged, the anisotropy of trans- 
verse conductivity also vanishes completely. This confirms 
that the combined effect of Rashba SOC, in-plane magne- 
tization, and nonmagnetic remote disorder could lead to 
AMR. Our numerical evaluation shows that these longitu- 
dinal and transverse AMRs are consistent with the stan- 
dard phenomenology due to symmetry arguments [6,18]: 
^<y xx /a av = C7cos2£, Aa yx /a av = C7sin2£. Here Cj is 
a dimensionless constant and is sometimes called noncrys- 
talline coefficient in literatures. One note that crystalline 
AMR coefficient vanishes in this case, which is a special 
property of the Rashba model. It is not valid for systems 
with other SOC, such as Dresselhaus SOC [18]. 

Now we limit ourselves to AMRs defined as the relative 
change between longitudinal resistivities for magnetization 
along and normal to the current direction. We take the 
current direction along [100]-axis. In this case, it is found 
that the transverse conductivity vanishes. Hence, AMR is 
given by 

AMR = — — ^2^ = 2-p — . (14) 

Here a xx and a xx are the corresponding longitudinal con- 
ductivities for M || J and M _L J with J as the cur- 
rent density. Also one find that this definition of AMR is 
disorder-density- independent . 

In Fig. EJa), we plot AMRs as functions of spin-orbit 
interaction constant for various magnetizations at fixed 
impurity distance s = 30 nm. It is seen that the magni- 
tude of magnetization can affect AMR strongly. A large 
AMR (~ 24%) can be observed for Rashba coupling pa- 
rameter up to 6 x 10 11 eVm at M = 2.5 meV. With the 
increment of Rashba SOC coefficient, AMR ascends and 
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may saturate at strong coupling. We also evaluate the 
dependencies of AMR on the magnitude of magnetization 
for various SOC constants in Fig. [3fb). With the rise of 
magnetization, AMR first increases, and then decreases. 
However, AMR is always positive. The strength of the 
spin-orbit interaction can affect both the value and the 
position of maximum AMR. For short-range nonmagnetic 
impurity, AMR vanishes completely. 

In order to investigate the density-related feature of 
AMR, in Fig. [3 AMRs are calculated for various electron 
density for fixed magnetization and disorder distance. It 
is clear that AMR is very sensitive to the electron con- 
centration. With the increasing density, AMR, arising 
from the electric remote scattering in SOC semiconduc- 
tor with in-plane magnetization, drops quickly. For ./V = 
5.0 x 10 11 cm~ 2 at large coupling constant, AMR ~ 0.45%. 
It is small but still measurable experimentally [6] . The in- 
set shows the dependence of AMR on electron density for 
a — 4 x 10 _11 cVm. This implies vanishing AMR in the 
limit of k + (8f.) m k-(8k) with k±(8k) as the two angle- 
dependent Fermi wave vectors. 

We now make some comments on the experiments to 
confirm the present results. Since we deals with AMR 
arising from nonmagnetic disorder, the nonmagnetic n- 
type InAs-based heterojunction can be well satisfied. The 
in-plane magnetization may be induced by an in-planc 
magnetic field. The magnitude of magnetic field corre- 
sponds to a magnetization M = 1 meV is 2.16 T (in InAs- 
based heterojunction, the effective g-factor g — 8 [27]). 
The Rashba SOC constant can be tuned by controlling 
the gate voltage. The usual Hall setup is well satisfied 
to measure this AMR. We note that Papadakis et al. re- 
ported the observation of AMR of two-dimensional holes 
in nonmagnetic GaAs [28]. Our present study may pro- 
vide a possible interpretation of that novel phenomenon. 
However, careful theoretical investigation should be made 
for p-type semiconductor systems. 

Conclusion. — In summary, AMR for two- 
dimensional electron systems with a Rashba-type 
spin-orbit splitting and an in-planc magnetization is in- 
vestigated for nonmagnetic impurity collision. It is found 
that the combined effect of SOC, in-plane magnetization, 
and electric remote disorder leads to AMR. The disorder 
distance can strongly affect the degree of anisotropy of 
longitudinal and transverse conductivity. The strong 
density-related character of AMR is also demonstrated. 
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